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Selection of state-insensitive or "magic" trapping conditions with ultracold atoms 
or molecules, where pairs of internal states experience identical trapping poten- 
tials, brings substantial benefits to precision measurements and quantum computing 
schemes. Working at such conditions could ensure that detrimental effects of in- 
evitable inhomogeneities across an ultracold sample are significantly reduced. How- 
ever, this aspect of confinement remains unexplored for ultracold polar molecules. 
Here, we present means to control the AC Stark shift of rotational states of ultracold 
diatomic polar molecules, when subjected to both trapping laser light and an exter- 
nal electric field. We show that both the strength and relative orientation of the 
two fields influence the trapping potential. In particular, we predict "magic electric 
field strengths" and a "magic angle" , where the Stark shift is independent of the DC 
external field for certain rotational states of the molecule. 

The advantage of using state-independent light traps for precision frequency measure- 
ments with ultracold atoms has been demonstrated in several experiments [TJ [2] . Applica- 
tions of this approach were analyzed in the context of optical atomic clocks and coherent 
control of atoms and photons within an optical cavity. For these applications the frequency 
of the laser beam that creates a far-off resonant optical dipole trap is chosen such that the 
AC Stark shift of the ground and one excited electronic atomic level are the same. In this 
way any optical atomic transition between these levels is unaffected by the trapping light. In 
Refs. [SHE] state-insensitive trapping conditions have also been found for microwave transi- 
tions in the atomic ground state by using a combination of the vector and tensor components 
of the AC Stark shift and an external magnetic field. 

Ultracold molecular systems possess unique proporties that are considered to make them 
potentially useful as tools for precision measurements [TJ [8] and quantum computing [9j [10] . 
Therefore it is desirable to extend the zero- differential AC Stark shift technique to these 
more complex systems. The idea of "magic" frequencies for vibrational Raman transitions 
in homo-nuclear Sr 2 molecules was first explored in Refs. [HJH2]. This molecule is proposed 
for a search of possible time variation of the electron-to-proton mass ratio. The AC Stark 
shift of a molecule is determined by the dynamic molecular polarizability a(u), which is a 
function of radiation frequency v and its polarization. 

Polar molecules have a permanent dipole moment and their levels can be shifted and 
mixed with one another by applying an external electric field. This opens up a new way 
to create "magic" trapping conditions for two rotational levels of the molecule. In the 
presence of an external electric field, J is not a good quantum number and all states, even 
the "rotationless" ground state, have an anisotropic polarizability [T3J . The anisotropy of 
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the dynamic polarizability of these levels manifests itself as a dependence on the relative 
orientation of the polarization of the trapping laser and the DC electric field. The combined 
action of these two fields can be a powerful tool to manipulate and control ultracold molecules 
trapped in an optical potential. 

The behavior described here has potential applications to several experiments that have 
been envisioned for diatomic polar molecules held in optical traps. The implications are 
particularly striking for the use of polar molecules in an optical lattice as quantum bits. 
As was described in Ref. [9], a pair of rotational states forms a suitable quantum bit. 
However, it was pointed out that this system is susceptible to decoherence due to intensity 
fluctuations in the optical trapping lasers, if the dynamic polarizabilities of these levels 
differ as is generally the case. This in turn leads to very stringent requirements on the 
laser intensity stability for such a system to be practical. As we will show, this limitation 
can be removed by adjusting the experimental parameters to guarantee that the dynamic 
polarizabilities of these states are equal. In particular, for this proposed system (where a 
spatial gradient of the electric field £ is required, so working at a magic electric field value 
is impossible) it should be possible to use light polarized at the "magic angle" relative to 
the static field in order to eliminate this potentially dangerous source of decoherence. 

The existence of magic electric field values is also of possible use for envisioned applica- 
tions where polar molecules in optical lattices are employed in novel types of many-body 
systems. This includes, for example, cases where the properties of long-range molecular 
interactions are tailored by a combination of a static electric field and resonant microwave 
fields coupling different rotational states [T4T - [T6] . Working at a magic electric field value in 
such systems could ensure that the inevitable inhomogeneities in the intensity of the trap- 
ping light across a large sample would not change the resonant condition for the microwave 
drive fields. Hence, working under such "magic" conditions might be necessary to implement 
proposals of this type. 

Motivated by these ideas, we calculate the near-infrared dynamic polarizability of various 
rotational levels of the v = vibrational state of the X X E + potential of the KRb and 
RbCs molecules, under the simultaneous influence of trapping electromagnetic and static 
electric fields. We calculate dynamic polarizability using the computational techniques 
developed in previous publications [TTHT9] . Our calculations with an external electric 
field are predominantly performed at a laser frequency of 9174 cm -1 (or wavelength of 
1090 nm), which corresponds to an often-used frequency to trap atoms and molecules in 
ultracold experiments. In the near infrared this laser frequency is sufficiently far away from 
molecular resonances of the excited electronic states that heating due to photon scattering 
is negligible. We focus on external electric field strengths up to 15 kV/cm, a value which 
should be experimentally accessible. 

We begin by studying the dynamic polarizability or AC Stark shift for the J = and 
J = 1 rotational levels of the v — X 1 S + state of KRb and RbCs in the near infrared, 
without an external electric field. The procedure used to determine the complex molecular 
polarizability has already been described in our previous paper [T7] . Figure [I] illustrates the 
results of this study in the absence of an electric field by showing the dynamic polarizability 
a as a function of laser frequency hv. The range of laser frequencies spans a technologically 
relevant near-infrared optical domain and includes the lasing frequency near 9174 cm -1 
used in Ref. [13]. The figure also shows that the dynamic polarizability of the molecular 
state depends on the rotational quantum number J and its projection M. The curves for 
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FIG. 1: Dynamic polarizability in the absence of an electric field and in atomic units for the v=0, 
J=0 and 1 levels of the X 1 E + ground state of KRb and RbCs as a function of trapping frequency 
in the near infrared domain. The dashed lines show the polarizability for the J=0 and M = 
level, which is independent of the light polarization. The solid lines correspond to the J= 1 state 
with magnetic sublevels M = and ± 1 illuminated by linear polarized light along the z direction. 
One atomic unit of polarizability corresponds to 4.68645xl0~ 8 MHz/(W/cm 2 ). 

the different states do not cross as a function of the laser frequency. In other words, at zero 
electric field we can not find a "magic" frequency for J=0 and J=l states using trapping 
light in the near infrared. 

Field-dependent dynamic polarizability 

Here we extend the idea of the AC Stark shift for the rotational levels of the X 1 E + ground 
state to the mixing of these levels in a static electric field £ = £z along the space-fixed z 
direction. Our analyses show that the experimentally accessible electric fields mix only a few 
low-lying rotational states. The ultimate goal is to introduce the dynamic polarizability of 
the mixed rotational levels as a function of an external DC electric field and the polarization 
of the AC trapping field. In a heteronuclear molecule this Stark mixing is primarily due 
to the permanent electronic dipole moment d; we neglect the much smaller effects due to 
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Stark mixing with other electronic states. For the Hund's ) X X S + state of an alkali- 

metal dimer, the total electron spin and orbital angular momentum are not coupled to the 
molecular rotation. The molecular wavefunction in the lab frame, \vJM) z , is then given by 

\vJM) z = \vJ) x Yj M (R) = |X 1 S + )| x Y JM (R) , (1) 

where v is a vibrational quantum number, J and M are the molecular rotational angular 
momentum and its projection along the z axis, ip v j(R) is the radial rovibrational wave- 
function, Yjm{R) is a spherical harmonic, R is the orientation of the molecule relative to 
the electric field direction z, |X 1 S + ) is the electronic wavefunction with projections defined 
along the internuclear axis. For the X 1 S + state, in even modest electric fields the nuclear 
spins are well-decoupled from the other spins and angular momenta; hence we ignore the 
nuclear spins here. Then the Hamiltonian for such a system becomes 

H=J2 E vj\vJM) z x z (vJM\ - d- £, (2) 

vJM 

where E v j = G v + B V J(J + 1) and G v and B v are the vibrational energy and rotational con- 
stant of vibrational level v, respectively. Higher order rotational corrections are negligible. 
For the alkali-metal dimers KRb and RbCs we have that AG V = G v+ \ — G v is on the order 
of 50-100 cm -1 for small v, while B v is on the order of 0.017-0.037 cm -1 . 

We evaluate the matrix elements of the operator — d ■ £ by noting that after av- 
eraging over the electronic wavefunction |X X S + ) it reduces to —d(R)C w (R)£ z , where 
d(R) = (X^+lfilX 1 !! 4 ") is the .R-dependent permanent electric dipole moment, Ci m (R) = 

^ An/ (21 + l)Yi m (R) are tensors of rank I, and £ z is the electric field strength. Consequently, 
this operator conserves the projection quantum number M. The matrix element between 
two rovibrational states is 

z (vJM\ - d ■ £\v'J'M') z = -6 M m> d vjyj'£z J dRYj M (R)Cw{R)Yj'M'{R) (3) 

= -5MM'd v j, v >j l E z (-l) M y/(2J + 1)(2 J' + 1) 

where 5mm 1 is the Kronecker delta function, d v j tV ij> = dRip v j(R)d(R)ilj v iji(R) 1 and (• • •) 
are 3-j symbols (see e.g. |20j). This matrix element is nonzero when J + 1 + J 1 is even, 
according to the parity selection rules, and is independent of the sign of M and M' . For the 
small J values of interest here we can assume that the J dependence of ip v j(R) is negligible. 
Moreover, coupling between vibrational levels can also be ignored, as the dipole moment 
d(R) is a slowly varying function with R and the spacing between vibrational levels is large 
compared to the rotational splitting. For tensor operators Ci m (R) of rank 0, 1, and 2 the 
relationship between Cartesian x, y, and z components and spherical m = —1, 0, and +1 
components can be found in Refs. . 

For each projection M and vibrational level v, eigenvalues and eigenvectors of the Hamil- 
tonian, Eq. ([2]), are obtained by the direct diagonalization of the Hamiltonian matrix in- 
cluding rotational states J from \M\ up to some value J max - For the external electric 
field strength and J values of interest, J max = 10 is sufficient for convergence. We la- 
bel the eigenenergies by E v j M with corresponding eigenvectors \v.JM) = U v j M j\vJM) z . 
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Here J is an integer index with values J = \M\,\M\ + 1,..., such that the eigenstate 
\vJM) adiabatically connects to the electric field-free eigenstate \vJM) z with J = J. For 
\M\ > the levels with projection quantum number — M and M remain degenerate. The 
dipole matrix element between states \vJM) and \v'J'M') of the X potential is given by 
(vJM\d a \v'J'M') = Ej,j>UfjU v j', A j', x z (vJM\d a \v'J'M') z , where a is a spatial index that 
can be expressed either in Cartesian coordinates x, y, z, or spherical coodinates q — 0, ±1. 

Now we are able to calculate the dynamic polarizability of the mixed rotational eigenstates 
at laser frequency v. It is determined by the properties of the operator a, ca i{y) defined by 

with 7 enumerating eigenstates of the ground as well as excited electronic potentials in 
the presence of the electric field and a, a' = x, y, or z. We will focus on the J = and 
J = 1 levels and without loss of generality assume that the x axis of our coordinate system 
lies in the plane spanned by the electric-field direction and the orientation of the linearly- 
polarized laser light. Moreover, we are interested in the situation where the level shifts due 
to the laser are small compared to those induced by the electric field. For the isolated levels 
with M=0, the polarizability is determined by the diagonal matrix element a v (T J J' I=Q {v) = 
(vJM = 0\a aa '(u)\vJM = 0). The \M\ > dynamic polarizability needs to be treated 
by degenerate perturbation theory within the two-dimensional subspace \vJM) and \vJ — 
M). In fact, our choice of the x direction and the symmetry properties of the dynamic 
polarizability ensure that the linear combinations \v JM, ±) = { \v JM) ± \v.J—M)}/\/2 with 

M > are the correct eigenstates. For these states the dynamic polarizability a^J*' (u) = 
(vJM, ±\a a(T i{v)\vJM, ±). For diatomic species, only the diagonal elements a^ M,± , a^ M,± , 

and a v z J z M,± are nonzero. Hence in an oscillating electric field £ (t) = £ o (0)Re{ee j27rv *} (where 
e*is the complex unit vector indicating the polarization), the state \vJM, ±) shifts in energy 

by an amount AE given by AE = —J2a,cr' l^o(0)| 2 a^i V/ ' ± (^)e CT e*,/4. The sum is over the 
spatial indices a, o' = x, y, z. 

We are interested in near infrared laser frequencies, which are detuned away from res- 
onances with rovibrational levels of the electronically excited potentials. In particular, we 
focus on wavelengths between 1000 nm and 1100 nm. Starting from the v — vibrational 
level of the X X S + states of KRb and RbCs, photons of this wavelength do not possess suffi- 
cient energy to reach rovibrational levels of the electronically excited singlet : A potentials. 
The photons do have enough energy to reach vibrational level of the triplet b 3 n potential 
considering only single-photon excitations. However, such transitions require relativistic 
spin-orbit coupling to the A 1 S + potential to acquire a non-zero dipole matrix element; such 
spin-orbit induced couplings are small enough to neglect under the conditions of interest 
here. 

This allows us to make several approximations. Firstly, we can use non-relativistic po- 
tentials and transition dipole moments and in the calculation of the dynamic polarizability 
only consider singlet 1 S + and x n potentials. Secondly, assuming a large detuning (such that 
\E 1 — E v j M — hu\ ^> E v jm for all states of interest) we can neglect the electric field and 
rotational dependence of the energy denominators in Eq. Q. Finally, we find that in the 
near infrared the contribution to the polarizability from intermediate states 7 in the ground 
X 1 S + state is small. 
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With this in mind the polarizability becomes 



vJM,± (i \ r^j jrv\M\ TT v\M\ 



(J(7 



^ J ' J J ' J> e tl I ^e, e A -E v +hv E eVeA -E v -hv\ K ' 

x ^ z ( U JM,±|4|e^e^M e A), z (e ? ; e J e M e A|4'k^M,±) z , 

J e M e 

where the energy E v is a typical vibrational energy in the ground state potential, and the 
rovibrational wavefunctions \ev e J e M e A) of the electronically excited states with approximate 
energy E eVe A are given by 

\ev e J e M e A) z = \ev e A) x \J e M e A) z = |M^| e *A) j x \^^Df* k {R) 

where the vibrational and electronic dependence has been isolated in the ket \ev e A) and 
the rotational dependence in | J e M e A) z , respectively. The wavefunction <j) Ve (R) is the radial 
rovibrational wavefunction, the ket |e X A) is the electronic state with projection quantum 
number A defined along the internuclear axis, and D J MM ,(R) is a Wigner rotation matrix 
that describes a symmetric top rotational wavefunction. The sum over A includes both 
positive and negative values, where A = corresponds to excited 1 S + electronic states and 
A = ±1 to X U states. 

As stated before the energy of the excited state and (fi Ve (R) depend on the electronic 
state \e l A) and vibrational level v e , but not J e and M e . Consequently, the transition dipole 
moments separate into z (vJM,±\d a \ev e J e M e A) z = (vJ\d(R)\ev e A)Fj M vI ^ A , where Fj^ A is 

an integral of the product of three Wigner rotation matrices D MM ,(R) over the orientation 
of the molecule R that can be evaluated using angular momentum algebra [20J. 

Moreover, we have verified that (v J\d(R)\ev e A) is nearly independent of J. The sums 
over J e and M e in Eq. ^ can now be performed and we finally find 

(vJM,±\a ca ,(y)\vJM,±) = (6) 

E™ A (,A TT v \M\ttv\M\ v pJeM e A rpJ e M e A 

A J,. J' ' J e M e 

with the A-dependent 

a\u) = J2(vJ\d\ev e A)(ev e A\d\vJ) {- 1 —— + ■ (7) 

The parallel oP[y) and perpendicular a l {v) contributions to the polarizability are due to 
transitions to the E and II states, respectively. 

Magic DC Electric Field 

The polarizability depends on the stength of an external electric field through the use 
of the unitary matrices Uj A j of Eq. (6). Figure |2j shows the dynamic polarizability of 
the ground states of the KRb and RbCs molecules as a function of the external electric 
field strength. The left panels of these figures correspond to the parallel (e* = z) and 
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FIG. 2: Dynamic polarizability at a wavenumber of 9174 cm -1 of the v=0, J=0 and 1 levels 
of the X^" 1 " ground state of KRb (top row) and RbCs (bottom row) as a function of the field 
strength of an external electric field. The polarization of the trapping field is parallel (left panel) 
and perpendicular (right panel) to the direction of the electric field. The circles indicate a crossing 
point between J=0 and 1 polarizabilities. The degeneracy of the states J = 1 and M = ±1 is lifted 
by perpendicularly polarized laser light. The new states are labeled | + 1) ± | — 1), corresponding 
to | J = 1,-M±) respectively. The M = ±1 states remain degenerate for parallel polarization. 



the right panels to the perpendicular (e* = x) polarization of the trapping light relative 
to the electic field direction 8. In all cases the polarizability depends on both the state 
index J and the angular momentum projection M. For KRb (Fig. [2, top row) a "magic" 
electric field strength exists at £ = 10 kV/cm, where the polarizability of J=0, M=0 and 
J—l, M=0 states coincide. This is possible due to the polar character of the molecules. 
For RbCs (Fig. [2j bottom rowj two "magic" electric field strengths exist. For the pair 
of states J — 0, M — and J — 1,M = a crossing occurs near 2 kV/cm. Another 
crossing appears for the states J = 1, M — and J — 1,M — ±1 at 4.7 kV/cm. They 
both are at a much smaller field strength than for KRb, since RbCs has a larger permanent 
dipole moment and smaller rotational splittings in the ground state. Note that the 
M = "magic" electric field occurs at the same field strength for both parallel and per- 
pendicular polarization of the trapping light. In fact, they are the same for any polarization. 
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FIG. 3: Dynamic polarizability of the v=0, J=0, M= level of the X 1 !]" 1 " ground state of RbCs 
as a function of external electric field strength and the angle 8 between the direction of the electric 
field 6 and the polarization e of the trapping light at a wavenumber of 9174 cm -1 . The contour 
lines are marked by the polarizability value in atomic units. 



Magic Angle 

In future experiments one would expect to be able to change the angle between the static 
and dynamic electric fields. Figure [3] shows a surface plot of the dynamic polarizability as a 
function of £ and the angle 9 between £ and polarization of the trapping light e, for linearly 
polarized light. The polarizability depends smoothly on both 9 and £ . 

Figurek compares the dynamic polarizability of the states | J — 0) and | J — 1, M — 0, ±1) 
as a function of the angle 9 of the linear polarization of the optical field relative to the static 
field (such that e = cos 9z + sin 9x), for several static electric field strengths within the range 
from to 6 kV/cm for KRb and to 3 kV/cm for RbCs. All curves with M = cross at 
the angle 9 = 9 such that cos 2 #o = 1/3? or 9o ~ 54 degrees. We refer to this as the "magic 
angle" since here the AC Stark shift is independent of the internal state of the molecule. 
This behavior occurs in many contexts and is a simple consequence of the rank-2 tensor 
structure of the polarizability [24J. The angular dependence of the J = level is smaller 
than that of the J = 1 levels as expected. It inherited the zero-electric field properties of 
the scalar J = state. As seen in Fig. [2J the RbCs polarizability of the J = 0, M = and 
J = 1,M = levels cross at a magic electric field strength of 2 kV/cm. For Fig. [4] this 
implies that the curves start to "overlap" when £ is equal or larger than this magic field 
value. 

The polarizability operator a (JO -/(z/), defined by Eq. Q, is a reducible rank-two tensor 
operator. Hence it can be expressed as a sum of irreducible tensor operators a^iy) of 
rank k = 0, 1,2. In terms of these irreducible tensor operators, the AC Stark shift AE is 
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FIG. 4: Dependence of the dynamic polarizability on the angle between the polarization e of the 
trapping light at a wavenumber of 9174 cm -1 and the direction of the external electric field £, for 
different values of £ . The polarizability of states with projections M = and M = ±1 of the X 1 !]" 1 " 
ground state of KRb and RbCs are shown in the left and right panels, respectively. For clarity 
£ ranges from 0.6 kV/cm to 6 kV/cm in steps of 0.6 kV/cm for KRb and from 0.3 kV/cm to 3 
kV/cm in steps of 0.3 kV/cm for RbCs. At the "magic angle" = 6q = 54°, the polarizability of 
the states J = 0, M = and J = 1, M = are identical and independent of the strength of the 
static electric field. For the J = 1, M = ±1 states no magic angle exists. 



proportional to the diagonal matrix element of the operator QwdV) 6 ^' , which can be 
written in the general form 

E E (8) 

k=0 q=-k 

where q is a spherical tensor projection index and the explicit forms of the irreducible 
spherical tensors a* (v) and ee^ are given in Refs [211 22] • 

In order to derive the "magic angle" condition for M = states, we consider the effect of 
each term in the expansion (|8]). The term with k = corresponds to the scalar polarizabil- 
ity; the operator a°(i/) = J2aa> oi a(T i{u)5 a a' has, under our approximations, diagonal matrix 
elements that are independent of J (or J) and M for all states of interest, and similarly the 
quantity eefj oc e • e* = 1 is independent of the polarization of the optical field. The term 
with k = 1, corresponding to the vector polarizability, in general is significant. However, 
for the special case of linearly polarized light where e*is real, the quantities ee^ = (e x e*) q 
vanish and hence the effect of the vector polarizability is zero. For the tensor polarizability 



10 



terms (with k = 2), from the Wigner-Eckhart theorem only the operator component with 
q = gives rise to a non-zero diagonal matrix element for M = states. Hence the contri- 
bution of this term is proportional to ee\. Without loss of generality we can define the linear 
polarization as e = e x x + e z z = cos 9z + sin 9x. In this case eel ^ e z e * — e - e**/3 = cos 2 9 — 1/3. 
Hence the contribution to AE due to the tensor polarizability also vanishes for all states 
with M = 0, when the optical field is linearly polarized at the "magic angle" 6 = 9q. Under 
this condition the only contribution to the dynamic polarizability is from the k = scalar 
term, which is the same for all states of interest. 

The property of the "magic" electric field discussed in Fig. [2] can also be understood 
in terms of the tensor structure of the polarizability. At certain values of the applied DC 
electric field E, the rank-2 components of the polarizability, a 2 ,, of the J = 0, M = and 
J = 1, M = levels becomes the same. When this condition is met, the AC Stark shift 
becomes independent of the direction of the trapping light's linear polarization. 

Here we finish with an example of a specific implementation of a "magic angle" 3-D 
lattice. Let the lattice be formed by three orthogonal retroreflected laser beams a, b, and c, 
with initial propagation directions k given by k a = y, kb = (x + z)/ \/2, and k c = (x — z)j \/2. 
These three beams each have a different frequency z/, such that v a = — 5b = v c — 5 C ; 
here the offset frequencies 5^ c must satisfy z/ a fe c 3> 5 fe c 3> f mot , where f mot is the motional 
frequency of the molecules in the optical trapping potential. The use of different frequencies 
(which can be generated from a single laser by using e.g. acousto-optic modulators) in 
this manner eliminates the effect of interference terms between the different laser beams: 
such terms average to zero rapidly over the time of motion of the atom, and hence can be 
neglected. The resulting average trap potential is then simply the sum of the potentials due 
to each individual laser beam. Finally, the polarizations of the three beams can be chosen 
as e a = ^2/3x + \Jl/3z; ib = j2/3k c + ^Jl/3y; and e c = ^2/3kb + v/l/3y. In each case, 
|e • z\ — cos 9 . 
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